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Abstract 

We discuss single diffractive gauge bosons (7*, W^, Z) production in proton-proton collisions 
at different (LHC and RHIC) energies within the color dipole approach. The calculations are 
performed for gauge bosons produced at forward rapidities. The diffractive cross section is predicted 
as function of fractional momentum and invariant mass of the lepton pair. We found a dramatic 
breakdown of the diffractive QCD factorisation caused by an interplay of hard and soft interactions. 
Data from the CDF experiment on diffractive production of W and Z are well explained in a 
parameter free way. 
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I. INTRODUCTION 



The characteristic feature of diffractive processes at high energies is presence of a large 
rapidity gap between the remnants of the beam and target. The general theoretical frame- 
work for such processes was formulated in the pioneering works of Glauber [l[ , Feinberg and 
Pomeranchuk 0] , Good and Walker [3| . The new diffractive state is produced only if differ- 
ent Fock components of the incoming plane wave, which are the eigenstates of interaction, 
interact differently with the target. As a results in a new combination of the Fock compo- 
nents, which can be projected to a new physical state (e.g. see review on QCD diffraction 
Ref. fl). 

The main difficulty in the formulation of a theoretical QCD-based framework for diffrac- 
tive scattering arises from the fact that it is essentially contaminated by soft, non- 
perturbative interactions. For example, diffractive deep- inelastic scattering (DIS), ^*p — > 
Xp, although it is a higher twist process, is dominated by soft interactions Within 
the dipole approach such a process looks like elastic scattering of qq dipoles of different 
sizes, and of higher Fock states containing more partons. Although formally the process 
7* X is an off-diagonal diffraction, it does not vanish in the limit of unitarity saturation, 
the so called black disc limit. This happens because the photon distribution functions and 
hadronic wave functions are not orthogonal. Such a principal difference between diffractive 
processes in DIS and hadronic collisions is one of the reasons for breakdown of diffractive 
QCD factorization based e.g. on the Ingelman-Schlein model 0]. In particular the cross 
section of diffractive production of the W boson was found in the CDF experiment (jl, 0] to 
be six times smaller than was predicted relying on factorization and HERA data |10j | . The 
phenomenological models based on assumptions of the diffractive factorisation, which are 



widely discussed in the literature (see e.g. Refs. [ll|, [12|), predict a significant increase of 
the ratio of the diffractive to inclusive gauge bosons production cross sections with energy. 
This is supposed to be tested soon at the LHC. 

The process under discussion, diffractive abelian radiation of electroweak gauge bosons, 
is the real off-diagonal diffraction. It vanishes in the black-disc limit, and may be strongly 
suppressed by the absorptive corrections even being far from the unitarity bound. The sup- 
pression caused by the absorptive corrections, also known as the survival probability of a 
large rapidity gap, is related to the initial and final state interactions. Usually the survival 
probability is introduced in the diffractive cross section in a probabilistic way and is esti- 
mated in oversimplified models, like eikonal, quasi-eikonal, two-channel approximations, etc. 
The advantage of the dipole approach is the possibility to calculate directly (although in a 
model-dependent way) the full diffractive amplitude, which contains all the absorption cor- 
rections, because it employs the phenomenological dipole cross section fitted to data. Below 
we explicitly single out from the diffractive amplitude the survival probability amplitude as 
a factor. 

Another source of factorization breaking is the simple observation that diffractive abelian 
radiation by a quark vanishes in the forward direction (zero momentum transfer to the 
target) Indeed, the Fock components of the quark with or without the abelian boson 
(7*, Z, W, Higgs boson) interact with the same total cross sections, because only the quark 
interacts strongly. Therefore, after integration of the amplitude over impact parameter the 
Fock state decomposition of the projectile remains unchanged, and only elastic qp scattering 
is possible. Notice that nonabelian radiation (gluons) does not expose this property, because 
the Fock components \q) and \q g), although have the same color, interact differently [l3, 14 
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In the case of pp collisions the directions of propagation of the proton and its quarks do 
not coincide. Already this is sufficient to get a nonvanishing diffractive abelian radiation in 
forward scattering. Moreover, interaction with the spectator partons opens new possibilities 
for diffractive scattering, namely the color exchange in interaction of one projectile parton, 
can be compensated (neutralized) by interaction of another projectile parton. It was found 



in Refs. [151 . |16| that this contribution leads to a dominant contribution to the diffractive 
abelian radiation in the forward direction. This mechanism, leading to a dramatic violation 
of diffractive QCD factorisation, is under consideration in the present paper. The breakdown 
of the diffractive (Ingelman-Schlein) QCD factorisation is a result of the interplay between 
the soft and hard interactions, which considerably affects the corresponding observables 
15| . Recently, such an effect has been analyzed in the diffractive Drell-Yan process (lol . 16]. 
and here we extend our study of this interesting phenomenon to a more general case - the 
diffractive gauge boson production. 



II. DIFFRACTIVE GAUGE BOSONS PRODUCTION AMPLITUDE 

Consider first the general formalism for diffractive radiation of the electroweak gauge 
bosons, 7*, Z°, W^, within the color dipole approach. Let us start with consideration of 
the distribution functions for the Fock states contributing to heavy gauge bosons radiation 
by a quark (valence or sea) in the projectile proton. 



A. Gauge bosons radiation by a quark 

The qj — > g/7* transition amplitude is given by the vector 7*g coupling only, i.e. 

T(q f ->• q f j*) = -ieZgE^Ufj^Uf , (2.1) 

where Zf is the quark charge, Uf and Uf are spinors for the quark of the flavour / in the 
final and initial states, respectively. 

The couplings of Z° and W bosons to quarks contain both vector and axial- vector parts. 
The qj — > qfZ° transition amplitude is given by 

T(q f q f Z°) = — — e» Xz u f [g Z vJ l» - g Z aJ lM^ (2.2) 
while the qf — > qfW ± amplitudes read, 

T (QU -»■ Qf d W + ) = -7= . V fufd e^Uf^jjf, - g aJ l^s]u u , (2.3) 



T (Qh -> <lfuW~) = -7= . V fdfu e" Xw u fd [g vJ ^ - g aJ i^]u fu , 
z V z sm v w 

for the up (f u = u, c, t) and down (f d = d, s, b) quarks, respectively. Here, Vf u f d is the 
CKM matrix element corresponding to f u — > fd transition, and 9\y is the Weinberg angle. 
The weak mixing parameter, sin 2 ^, is related at the tree level to Gp, M z and a em by 
sin 2 9w = 47ra em /v / 2Gi?M| (we adopt here a em (m z ) — 1/127.934). The vector couplings 
at the tree level are 

9lu = \-\^ 2 6 Wl ^ /d = _I + ^ si n 2 ^, g% f = l] (2.4) 
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whereas axial- vector couplings are 



<U = \ < fd = -\, 9^ = 1- (2.5) 

Heavy gauge boson polarization vectors describing transverse (T), Ac = ±1, and longi- 
tudinal (L), X G = 0, polarization states are defined in light-cone coordinates 1 as 

ex=± = (0, 0, e A=± ), e x=± = T-^(l, ±i), (2.6) 

£ « = [h °J ■ < 27 > 

Note that we work in the physical (unitary) gauge. The calculations are performed in the 
high energy limit, i.e. in the limit where q + is much larger than all other scales. 

Let us start with the radiation of a heavy gauge boson by a quark interacting with a proton 
target. We assume that the longitudinal momentum of the projectile is not changed signifi- 
cantly by the soft interaction at high energies. In the high energy limit the corresponding s 
and w-channel amplitudes of the gauge boson bremsstrahlung in the quark-target scattering 



can be written as follows (cf. diffractive DY amplitude in Ref. |16| ) 
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where G = 7, H^ 1 * 1 , Z is the gauge boson under consideration; r? 2 = (1 — a)M 2 + a 2 m 2 ; 
a is the fractional light-cone momentum carried by the gauge boson, which has invariant 
mass M. The vector (v) and axial-vector (a) couplings g v j a j are defined in Eqs. (12 .4p and 
(12. 5p . whereas our notations imply , = 1 and g^/ = 0> °i,2 are the helicities of initial 
and final quarks, respectively; k± = p%± — pu_ + q± is the transverse momentum of the 
exchanged gluon; and l± = f>2± — (1 — cx)q±/ a is the transverse momentum of the final quark 
in the frame, where z-axis is parallel to the gauge boson momentum. The amplitude A for 
scattering of the quark on the nucleon in the target rest frame has the following approximate 



form, [23 



A aai (k ± ) ~ 2p u 1 5 aai V q (k 1 _), A a(J2 {k x ) ~ 2p J 2 5 tTtT2 V q (k ± ) , 

where the factorized universal amplitude V q (k±) does not depend on energy and helicity 
state of the quark. The coupling factor Cf, f = q,l, introduced in Eq. (12.81) . is defined for 
G = 7*, Z° and W ± bosons, respectively, as 



✓77 _ rz v nZ _ \/ Q ™ pW+ _ \/ a em T/ n w~ _ \/ a em T/ 

c f -v^z f , c f -—^-, Cj - 2 ^^ g v fufd , c f - 2 ^ 6 Jm^ 



1 As usual, the light-cone 4- vector p is defined as p — (p + , p , p), where p^ = p° ± p z 
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where a em = e 2 /(47r) = 1/137 is the electromagnetic coupling constant. In the case of gauge 
boson couplings to leptons we should substitute Vf u f d = Vf d f u = 1. 

Eventually, we can switch to impact parameter space performing the Fourier transforma- 
tion over l± and k±, and write down the total amplitude M q for gauge boson radiation in 
quark-proton scattering as follows, 



M£(b, r) = -2ipl K- A (f, a, M) ■ V q (b) - V q (b + ar) 



a 
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V q (b) = J ^ e -^V# ± ), ^(r, a, M) = V z v > w (r, a, M) - *% w (r, a, M) , 

where ar is the transverse separation between the initial and final quarks; tyy, A (f, a, M) 
are the light-cone distribution functions of the vector q — > Vq and axial-vector q — > Aq 
transitions in the mixed representation defined as, 

K(r a M) — C G q G a 3 VT^ I f!^ e ~^ ^ 2 ( P f)rMP2 + q) 

mr aM)= C G q G a 3 VT^ [ f!^ e -^^Ml^M + g) 
•* A yi,u,M) SSm" vi u j ^n) 2 a 2 l 2 ± + r) 2 

For an unpolarized initial quark the interference terms between the vector and axial-vector 
wave functions cancel each other, i.e. 

E *v-a(<x,Pl)*v-a(<x,P2) = ^(a,pi)^(a,p 2 ) + ^(a,p!)^(a,p 2 ) . (2.9) 

Cl,<T2 

The bilinear combinations of the vector V and axial-vector A light-cone distribution func- 
tions, corresponding to radiation of longitudinally (A = 0) and transversely (A = ±1) 
polarized gauge bosons have the form 2 , 

4$(a, A) = E \ E <(A)^(«,Pi)e,(A)^*(a,p 2 ) (2.10) 

A=±l cri,(T2 

= C ^f^W q a'K ( VPl ) K ( VP2 ) + [l + (1 - a) 2 ] if^Kx ( V Pi) Ki ( W2 ) 1, 

I Pip 2 I 

*^(a, A)*?(a, A) = ~ E <( A = 0)^(a,p' 1 ) e ,(A = 0)¥?(a, A) 
" S^M 2 (1 - a) 2 Ko (^x) K ( W2 ) . 



7T^ 

^(a,p\)^>,p 2 ) = (2.11) 

P2( n G \2 ( - 
_^q\ya,q) I 2 2/ n „ \2 T ^ L„ h , /i „ „2 Pi ' P2 1 



, (2 - a) 2 K (VPi) Ko (r?p 2 ) + [l + (1 - a) 2 ] rf^-^K x ( VPl ) K, (r]p 2 

Z7T I pip 2 

* L A {oc,p x )^ L A \oc,p 2 ) = S^^jr/ 2 K (r / p 1 )K (r ? P2) + « 2 m;^K 1 ( VPl ) K x (^p 2 ) 
where the averaging over helicity of the initial quark is performed. 



In the case of a heavy photon 7* bremsstrahlung by a quark such formulae were derived in Refs. 221424 1. 
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B. Forward diffractive radiation from a dipole 



The amplitude of diffractive gauge boson radiation by a quark- antiquark dipole does not 



vanish in forward direction, unlike the radiation by a single quark [13] , |15| . This can be 
understood as follows. When a quark fluctuates into a state \qG) containing the gauge 
boson G, with the transverse quark-boson separation r, the quark gets a transverse shift 
Ar = af. Nevertheless, both Fock states \q) and \qG), interact with the target with the same 
total cross section, this is why G is not produced diffractively in the forward direction. The 
situation changes if the boson is radiated diffractively by a dipole. Then the quark dipoles 
with or without a gauge boson have different sizes and interact with the target differently. 
So the amplitude of the diffractive gauge boson radiation from the qq dipole is proportional 
to the difference between elastic amplitudes of the two Fock components, \qq) and \qqG) flo|. 
i.e. 



\ [,,„ ( r„. r. n ) = - 2zp° ^ ^ 2 " r) 2L . i f, , ( h. /• , ) - 2In. /, ,(/,./■„ + ,» F) (2.12) 



where F p is the transverse separation of the qq dipole. The partial elastic dipole-proton 
amplitude is normalized to the dipole cross section, which is parameterized by the following 



simple ansatz 32 



J d%2lmf el (b,f p ) = a (l -e- r 2/W) (2.13) 



This saturated form, although is oversimplified (compare with (33j), is rather successful 
in description of experimental HERA data with a reasonable accuracy. We rely on this 
parametrization in what follows, and the parameters <To and Rq(x), as well as the explicit 
form of the amplitude f e i(b, r), will be specified later. 

The diffractive amplitude f)2.12p . thus, occurs to be sensitive to the large transverse 
separations between the projectile quarks in the incoming proton. These distances are 
controlled by a nonperturbative scale, which is one of the reasons for the breakdown of 
diffractive QCD factorisation in the diffractive gauge bosons production (for more details, 
see Refs. [lEGi})- 



III. SINGLE DIFFRACTIVE CROSS SECTION 



The differential cross section for the single diffractive di-lepton (11 pair in the case of 7*, Z 
and Ivi pair in the case of W^) production in the target rest frame can be written in terms 
of the gauge boson production cross section at a given invariant mass of the di-lepton M. 
Integrating the cross section over the solid angle of the lepton pair and the boson transverse 
momentum q± we get for the diffractive Drell-Yan cross section [16], 



d 6 VL,T(pp -> pUX) 



a. 



d 5 a L , T (pp ->■ jry*X) 



(Pq^dxx dM 2 d 2 5 ± 3vrM 2 d 2 q ± dx! d?8 A 



(3.1) 



where x\ is the fractional light-cone momentum of the di-lepton, 5± is the transverse mo- 
mentum of the recoil proton, and q± is the transverse momentum of the outgoing photon (or 
di-lepton). Compared to our previous study, here we are going to look at the g^-dependence 
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of the diffractive DY cross section in a much wider range of di-lepton invariant masses 
accessible at the LHC 

In the case of the diffractive production of G = Z°, W bosons, it is convenient to employ 
the simple and phenomenologically successful model for the invariant mass distribution in 
the decay of an unstable particle (for details, see e.g. Refs. [2l|) and to present the differential 
cross section in the factorized form, 

1 Wl)X) = Br ' G -» »°M ^IT X) ■ (3-2) 

a z q±axi dM A d z o± d z q±dx\ d z o± 



where Br(Z° ->- £ l=w ll ) - °- 101 and Br(W ± J2i= e ,»,r - °- 326 are the leptonic 
branching ratios of Z° and W ± bosons, and po{M) is the invariant mass distribution of the 
dileptons from the decay of the gauge boson G, 

1 MTr(M) 
** M) = * (KP- m lAMT G (M )r ^M)IM«1. (3.3) 

Here, vrtQ is the fixed on-shell boson mass and Tg(M) is its total decay width defined in the 
standard way by substitution m G — > M, i.e. 



^ ,, rS 3a em M _, .„,. a Pm M 

T W (M) ~ - , em 2n , T Z (M) 



4 sin 2 6 W y ' 6 sin 2 26 



w 



160 



sin 9 W - 40 sin 2 9 W + 21 



(3.4) 



Let us assume that the gauge boson is emitted by the quark q\. As a result of the hard 
emission the quark position in the impact parameters, being initially f\, gets shifted to 
r*i + af. Applying the completeness relation to the wave function of the proton remnant in 
the final state 

+ af, f 2 , f 3 ; {x l q 2 -}, {xj' 2 -})^(f [ + af, f 2 , f 3 ; {x' 1 /-}, {x']' 2 '-}) 

f 

= sfr - f[ + a(f- f'))S(f 2 - f' 2 )5(f 3 - f 3 ) J] 5(x j q/g - x'{ /g ), (3.5) 

j 

where f, x l q i g are the transverse coordinates and fractional light-cone momenta of the va- 
lence/sea quarks and gluons, we get the diffractive G* production cross section in the fol- 
lowing form [l5|, ltj 



d 5 a\ G {pp ->• pG*X) 



i-r — V / d 2 r x d 2 r 2 d 2 r 3 d 2 rd 2 r' d 2 bd 2 b' dx q TT dx l dx\ 
nY QA-k 2 x\ J - L .- L y y 

' n — 11/1/ oct i 



d?qx_dx\ d 2 5± (2tt . 

q=val, sea 

x V X v G _ A (f, a, M)^ G _* A (f, a, M) \%[f u ? 2 , r 3 ; x q , {xf-}, {xf-})\ 2 

x A(fx, f 2 , f 3 ; b; f, a)A{f,f 2 , f 3 ; b'; f, a) e iS± <^ e ^< r -*"> (3.6) 

where is the proton wave function, the summation is performed over all valence/sea 
quarks and gluons in the proton, and the light-cone fraction of the quark emitting the gauge 
boson x l q = x q is fixed by the external phase space variables X\ and a due to the momentum 
conservation, namely, 

Xq = it ' Xl = P+ (3 ' 7) 
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where P\ is the 4-momentum of the projectile proton, q is the 4-momentum of the produced 
gauge boson, and 



— * — * 

A = -21m fei(b, n - r 2 ) + 21m f a (b, f i - r 2 + ar) 
-21m / eZ (6, fi - r 3 ) + 21m / e /(6, f i - f 3 + ar) 



(3.8) 



is the properly normalized diffractive amplitude, where f e i(b,fi — f 2 ) is the partial elastic 
amplitude for dipole of transverse size r colliding with a proton at impact parameter b 
to be specified below. As expected, the diffractive amplitude A is proportional to the 
difference between elastic amplitudes for the dipoles of slightly different sizes. This difference 
is suppressed by absorptive corrections, the effect sometimes called survival probability of 
large rapidity gaps. 

The amplitude Eq. (13.81) is the full expression, which includes by default the effect of 



absorption and does not need any extra survival probability factor 3 . This can be illustrated 
on a simple example of elastic dipole scattering off a potential. The dipole elastic amplitude 
has the eikonal form, 



Im f e i(b, r 1 -f 2 ) = l- exp [ix( 



n, 



%X{r 2 



where 



X(b) 



dz V(b, z), 



(3.9) 



(3.10) 



and V(b, z) is the potential, which depends on the impact parameter and longitudinal coor- 
dinate, and is nearly imaginary at high energies. The difference between elastic amplitudes 
with a shifted quark position, which enters the diffractive amplitude, reads, 

Im f d (b, f x - f 2 + ar) - Im f d (b, ft - f 2 ) = exp [z'x(ri) - ix(r 2 )] exp [iaf ■ Vx{n)] ■ (3.11) 

The first factor exp[ix(r*i) — ix(^*2)] is exactly the survival probability amplitude, which 
vanishes in the black disc limit, as it should be. This proves that the cross section Eq. (13. 6p 
includes the effect of absorption. Notice that usually the survival probability factor is 
introduced into the diffractive cross section probabilistically, while in Eq. (13.61) it is treated 
quantum-mechanically, at the amplitude level. 

All the elastic amplitudes in Eq. ( 13. 8p implicitly depend on energy. They cannot be 
calculated reliably, but but are known from phenomenology. Since large dipole sizes 
\fi — fj\ ~ b ~ R p , i j (R p is the mean proton size) are important in Eq. (I3.8p . the 
Bjorken variable x is ill defined, and the collisions energy is a more appropriate variable. A 
parametrization of the dipole cross section as function of s was proposed and fitted to data 
in Ref. [l^], and the corresponding partial dipole amplitude is given by (25-27 



lraf e i(b,f p ,s,x q 



— 2 exp 



00 is 




exp 



[b+fp(l-Xg)f 



Ri(s) 



8irB(s] 

[b + f p (l/2-x q )} 2 
2B(s) 




+ exp 



[b + f p x q 
2B(s) 



R 2 N (s) + R 2 (s)/8 



(3.12) 



3 Such a statement has already been made in a similar analysis of the diffractive heavy flavor production 
performed in Ref. 14[ and in our previous work on diffractive DY study 16 1. 



s 



where x q is the quark longitudinal quark fraction in the dipole defined in Eq. (13. 7p . and 
R (s) = 0.88fm( So /,) - 14 , R%(s) = Bj(s) - - 1 -(t 2 c X , 

Ms) = aZ(s)(l + ^4). (3.13) 



Here, the pion-proton total cross section is parameterized as 29[ cr^(s) = 23.6(s/s ) ' 08 mb, 
s = 1000 GeV 2 , the mean pion radius squared is (3l[ (r 2 ^ = 0.44 fm 2 , and the Regge 
parametrization of the elastic slope B^f(s) = B + 1ol f ln(s//i 2 ), with B = 6 GeV -2 , a'p = 
0.25 GeV" 2 , and /i 2 = 1 GeV 2 can be used. We employ the s-dependent parametrization 
(I3.12p in what follows, because diffraction is essentially controlled by soft interactions. 

Finally, we parameterize the proton wave function assuming the symmetric Gaussian 
shape for the spacial valence quark distributions in the proton, as 

\%(r 1: f 2 ,f 3 ;x q , {xf '-)| 2 = 3 -^e-^^p(x q , {*f -}) 

x 5(r 1 + f 2 + r 3 )S(l-x q -^2x j q/g ), (3.14) 

j 

where sum is taken over all valence/sea quarks and gluons not participating in the hard 
interaction, x q is defined in Eq. (13.71) . a = (r 2 h ) _1 is the inverse proton mean charge radius 
squared; p is the valence quark distribution function in the proton. Notice that this dis- 
tribution has a low scale, so the valence quark carry the whole momentum of the proton, 
while gluons and the sea are included in the constituent valence quarks. The Gottfried sum 
rule based on this assumption is know to be broken [30J, but we neglect the related ~ 20% 
correction. 

Integrating over the fractional momenta of all partons not participating in the hard 
interaction we arrive at the single valence quark distribution in the proton, probed by the 
hard process - radiation of a heavy gauge boson, 

j dxldxl 5(1 -x q -J2 K' 3 '-}, {s? 3 -}) = PM » (3-15) 

where q denotes the quark flavor emitting the gauge boson G with the fraction x q given by 
Eq. (13.71) . In the case of diffractive Drell-Yan reaction [l6|, generalization of the three-body 
proton wave function (13. 14j) including different quark and antiquark flavors leads to the 
proton structure function as, 

J2 Z2 M X *) + PM] = ~ F ^) ■ (3-16) 

Xq 

However, in the case of diffractive W and Z production the coupling factor C q g^, a varies 
for different (valence/sea) quark species in the proton, so one has to deal with the original 
quark densities. Similar to the diffractive DY case, in actual numerical calculations below, 
when summing up the contributions of different quark flavors, we will generalize the above 
approach including the sea quark and antiquark densities in the proton at the hard scale 
imposed by the mass of the gauge boson. Also, the interference terms between amplitudes 
corresponding to gauge boson radiated by different valence quarks separated by large trans- 
verse distances in the proton are strongly suppressed in the hard limit r <C Ro{s), and are 
neglected. 
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IV. SINGLE DIFFRACTIVE CROSS SECTION IN THE FORWARD LIMIT 



Since the fluctuation of the dipole size caused by the hard vector boson production, ar ~ 
a/(l — a)M, is typically much smaller than the (semi)soft hadron scale \f\ — fj\ ~ b ~ R p , 
i 7^ j, one can derive an approximate analytical formulae for the diffractive cross section 

Im f el (b, 4- + ar) - Im /„(&, R l3 ) ~ d J^l^lM ^ (4 .1) 

where R\j — fi — fj. For the sake of convenience, let us transform the integrals in Eq. (13. 6p 
by introducing new variables r 2 — > Rn and r 3 — > Ri 3 , so that, 

J d 2 ri d 2 r 2 d 2 r 3 e-^+r^Sin + f 2 + r 3 ) = ^ / d 2 i? 12 d 2 i? 13 e-^ + ^3+Mi 3 ). ( 4 . 2 ) 

Since in the forward limit 8± — > the 6-dependence comes only into the partial dipole 
amplitude f e i defined in Eq. (I3.12p . it can be easily integrated fl5| . 



d% gIm ^fr4) = ^ e -^/^ (S ) ; (43) 

with the energy dependent parameters defined after Eq. (I3.12p . 

We see that the amplitude of diffractive gauge boson emission in the dipole-target scat- 
tering (I2.12p integrated over b, 

[ d 2 bM qq (b,R tJ ,f,a) oc«^4(f -R^e-xy^, (4.4) 

is proportional to the product of the hard scale r ~ 1/(1 — a)M and the soft hadronic scale 
Rij ~ Rq ~ 1/Aqcd- This means that the single diffractive cross section depends on the 
hard scale as a S( i ~ r 2 ~ 1/M 2 . 

It is well-known that the cross section of diffractive Deep-Inelastic-Scattering (DDIS) 
oddis ~ t 4 is dominated essentially by soft fluctuations at large r (for more details, see e.g. 
Ref. Q), as correctly predicted by the diffractive (Ingelman-Schlein) QCD factorisation. 
This happens since the end-point qq dipole fluctuations, driving the cross section at a — > 
or 1, have no hard scale dependence for light quarks m q Q 2 . In this case, the Q 2 - 
dependence comes only into their weight as ~ 1/Q 2 , even though it is of the higher twist 
nature. 

In opposite, the single diffractive gauge bosons production cross section behaves as ~ f-R, 
soft and hard fluctuations contribute in this process on the same footing, and their interplay 
does not depend on the hard scale, similar to the inclusive gauge bosons production. Hence, 
the forward diffractive abelian radiation turns out to be of the leading twist nature, and 
the diffractive-to-inclusive production cross sections ratio can depend on the hard scale only 
weakly through the x-dependence of the saturation scale, or more precisely ^0(^2), and can 
only increase (see below). 

However, if one uses the conventional diffractive factorisation scheme 0] the single diffrac- 
tive cross section, similarly to the DDIS process, one does not find any soft-hard interplay 
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as observed above, and the cross section turns out to behave as ~ r 4 , providing the higher 
twist nature of the single diffractive process. Correspondingly, this strongly affects the M 2 - 
dependence of the diffractive-to-inclusive boson production cross sections ratio, such that it 
decreases with M 2 , as opposite to our observation above. 

Therefore, the fundamental interplay between the hard and soft interactions in the for- 
ward diffractive Abelian radiation is the major reason for the diffractive QCD factorisation 
breaking leading to quite unusual features of the corresponding observables (for a similar 
discussion in the diffractive DY, see Refs. [IE IS)- 

As we have emphasized above, this 
interplay is absent in the DDIS and in diffractive QCD factorisation-based approaches to 
the diffractive DY (see e.g. Ref. [HI ]) leading to the energy and scale dependence of the 
corresponding cross section which is completely opposite to the one predicted above by the 
Color Dipole model. 

Further, the integrations over i? 12 and can be performed analytically leading to the 
diffractive cross section (13.61) in the forward limit 5± — > 0, 



dV Ao (pp->j>G*X) 



where 



d 2 q±dxi d5 
i 



=o 24vr 3 R*{s) A 2 



da 



Pi 



2 / 



2a 



+ 



R 2 o( 



d rd r 



[r ■ r 



2a 
~3~ 



{A 2 -AA 1 f 



+ 



A\ 



(A 2 



AA 2 ) 2 



(4.5) 



*F A (f,a,M)* 



v-a\ 



f ', a, M) e 



iq ± -(r-r') 



A, 



2a 

Y 



Rn 



(4.6) 



Assuming gaussian <5j_ dependence of the cross section, the 5± integrated and forward 
cross sections are related as, 



da(pp -)■ pG*X) 1 d 3 a(pp -» pG*X) 



d 2 q±dxi 



B sd (s) 



d 2 q±dx\ d5 2 



(4.7) 



ch 



)/3 + 2ot'jp ln(s/so), is similar 



The slope of the single- diffractive cross section, B s d(s) t 
to the one measured in diffractive DIS. In the next section we will explicitly derive the 
diffractive slope from the explicit parameterisation for the partial dipole amplitude (I3.12p . 

Finally, one can explicitly calculate the remaining integrations in the transverse plane 
over f and r ' by means of the following Fourier transforms 



d 2 rd 2 r' 



K (rjr)K (r]r') e 



iq±-(r-f') 



16?H 



(v 2 + q\Y 

4 i 4 

rf + q\ 



d rd r 



, (r • r ') 



f\2 



Ki(r]r)Ki(r]r') e 



iq ± -(r-r') 



57T 



rf{rf + qjj 4 

And we arrive at the following expressions for the cross section of transversely and longitu- 
dinally polarized gauge boson production, respectively, 



d 4 a T {pp^pG*X) 



d 2 q±_dx\ 



+ 



A\ 



B s d{s) 24-7T 3 R%{s) A 2 i(A 2 — 4Ai) 2 (A 2 -AA 2 ) 2 



(4.9) 



2tt 2 



da 



X'l 



a 



Pq[— +Pq[ — 



a 



9 2 

m q a 



(^JV + (C) 2 (2 - a) 



G \2, 



Jl + 



+ [l + a-c) 2 ]^}; 



11 



d 4 a L (pp->pG*X) 



a 2 ctq(s) 1 



d?q±_dx\ 
(C?) a ^ 



+ 



^2 



B sd (s)2A^Rt(s) A 2 l(A 2 -AA 1 ) 2 {A\ - 4A|) 2 



x (4.10) 



7T 



M — J +p? — 



•'•1 



a: 



G4) 2 M 2 (1 - a) 2 + G4) 2 ^, 



{9a,q) a m qJp J 2f 



These expressions for the differential distributions in the transverse momentum of the pro- 
duced gauge bosons allow us to perform gl-integration via the substitution, 



Ji(q±,V) -> h(v) = J d 2 q ± Ji(q ± ,r]) 
Mq±,v) -> h(v) = J d 2 q ± J 2 (q±,r]) 



3r/ 4 ' 
16vr 3 
~3rj r 



The left integrations over a and x\ can be done numerically. 



V. DIFFRACTIVE VS INCLUSIVE PRODUCTION OF GAUGE BOSONS 



The dipole description of inclusive Gauge boson production can be obtained generalizing 
what is known for the inclusive Drell-Yan process 22|, |24|, |28[ . The cross section of inclusive 
production of a virtual gauge boson G* with mass M and transverse momentum q± has the 
form, 



d 4 a Xa ( PP ^G*X) 1 



d 2 q± dx\ 



Ef 1 da 



■t-1 



Pq[ — j +Pq( — 

x a J V a 



x 



(5.1) 



d 2 rd 2 r' -\a{ar) + a(ar') - a{a\r- r '|)| ^ G _ A (r, a, M)^ G _* A (f ', a, M) e ^<^ . 

The principal difference of the inclusive gauge boson production from the diffractive one 
is in the typical size of the dipoles involved in the scattering. As is seen from e.g. Eqs. ( 12.13!) . 
( 13. 8j) . the diffractive scattering is dominated by large dipoles scattering at the hadronic scale, 
with the transverse size r p = Rij ~ Rq (soft scattering), whereas the inclusive production 
cross section (15. ip is totally driven by small-size dipoles scattering with r p = ar <C R (hard 
scattering). Therefore, different parameterizations for the dipole cross sections must be 
used - in the diffractive case above we have adopted the KST parametrization for the dipole 
cross section (or the partial amplitude (13 . 1 2[) ) with s-dependent parameters introduced in 
Eq. (I3.13P (l3l 13], whereas in the inclusive production case the Bjorken x-dependent GBW 
parametrization Eq. (I2.13P [32| is better justified: 



cr = 23.03 mb, R = Rq{x 2 ) = 0.4 fm x (x 2 /x ) 



0.144 



£0 



3.04 xlO" 4 , (5.2) 



where x 2 = q /P 2 , with P 2 being the 4-momentum of the target proton. 
In the leading regime of ar, ar' <C Rq. 



-<{ a(ar) + a(ar') — cr(a\r — f'\) 



a 2 a 



r • r 



(5.3) 
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so the inclusive gauge boson production cross section at forward rapidities {x\ ^> x 2 ) reads, 



d 4 a Xo (pp^G*X) 



E 



da 



X'l 



d 2 q ± d Xl (2tt) 2 R%(x 2 ) 
d 2 rd 2 r' (r • r ') x (r, a > M)^ G _* A (f ', a, M) e ^-<^ . 



x 



(5.4) 



We observe that the integrals over a and r, r' have the same form as in the diffractive cross 
section Eq. (14. 5p . 

The M-dependence of the differential cross sections for di-lepton inclusive production via 
an intermediate photon 7* or a gauge boson G* can be presented similar to the diffractive 
case, as 12811. 



da 



Q<- 



da(pp -> 7*X) 



d?q±dxi dM 2 
da(pp (G* II, lvi)X) 
d 2 q±dxi dM 2 



3irM 2 
Bv(G^ll 7 ln)pa{M) 



d 2 q±dxi 

da(pp G*X) 
d 2 q^dx\ 



(5.5) 



where the resonance mass distribution pc{M) is given by Eq. (13. 3ft . 

Eventually, we arrive at a simple form for the ratio of the diffractive and inclusive cross 
sections for di-lepton production, 



da{ d c /d 2 q 1 _ dx x dM 2 

And 



a 2 R 2 {M 2 Jx lS ) a 2 ( 



da^jd 2 q 1 _dx x dM 2 6tt B sd (s) a 



1 



(A 2 - AA,) 2 



+ 



A 2 



(A 2 



AA 2 ) 2 



(5.6) 



where functions ^1,2,3 were defined in Eq. (14. 6p . and fraction x 2 is explicitly given in terms 
of other kinematic variables in Eq. (14.61) . 

It turns out that the ratio (15.61) does not depend either on the type of the intermediate 
boson, or on its helicity Ac- To a good approximation, it is controlled mainly by soft inter- 
action dynamics, in terms of the soft parameters only Rq, Rq, a and a . A slow dependence 
of these parameters on the collision energy s, the hard scale M 2 and the boson transverse 
momentum q± completely determines such dependence of the diffractive-to-inclusive pro- 
duction ratio. A measurement of the M 2 (or q±) dependence of this ratio would allow to 
probe the x-evolution of the saturation scale, as well as to constrain its energy dependence. 
Hence, such a quantity is a very useful probe for the underlined QCD diffractive mecha- 
nism and the saturation phenomenon, and will be quantified based on existing KST/GBW 
parameterizations in the next section. 



VI. BREAKDOWN OF DIFFRACTIVE FACTORIZATION 

It is instructive to trace the origin of QCD factorisation in inclusive processes within 
the dipole description. The 1/Q 2 dependence of the DIS cross section at small x originates 
from two different sources. Most of the qq fluctuations of a virtual photon have a small 
size, r 2 ~ 1/Q 2 , except the endpoint (aligned jet) configurations with a — > 0, 1. The latter 
have large hadronic size and cross section, but their weight is small ~ 1/Q 2 - Thus, both 
contributions to the cross section behave as 1/Q 2 - 
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Similarly, in the Drell-Yan process of radiation of the heavy photon with fractional mo- 
mentum x\ the mean size of the artificial dipole [22| has small size r ~ 1/(1 — az)Q, except 
the endpoint configurations with a — > 1. Like in DIS, a is not an observable, but is inte- 
grated from x\ to 1 (see Eq. (15.1 1) ). This similarity, which reflects the factorisation relation 
between the two processes, also demonstrates its limitations. At large X\ — > 1 the inclusive 
Drell-Yan reaction is fully dominated by the soft component and factorisation breaks down. 

For the diffractive channels the factorisation relation breaks down at any x\. DIS diffrac- 
tion is fully dominated by the soft dynamics, since the probability of endpoint configurations 
is still the same, oc 1/Q 2 , while the cross section is enhanced by a factor of Q A compared 
with the hard component 0] . However, the mechanism of Drell-Yan diffraction is quite dif- 



ferent [15j, [16|, apparently breaking the factorisation relation. It comes from the hard-soft 
interference, which imitates a leading twist throughout the whole range of X\. Again, like in 
inclusive process, the hard and soft (endpoint) components make comparable contributions 
to the diffractive Drell-Yan cross sections. Of course, this is true for other gauge bosons as 
well. 

Although involvement of large distances in diffractive heavy boson production is in obvi- 
ous contradiction with factorisation of hard and soft scales, an observable manifestation of 
that is not trivial. Indeed, the cross section of a hard process qq — > // with the sea quark 
density in the Pomeron measured in diffractive DIS, although it is a higher twist, imitates 
the leading twist scale dependence. Nevertheless, the predicted scale and energy depen- 
dences are quite different, as we demonstrate in the next Section IVIH Notice that a much 
more pronounced breakdown of diffractive factorisation was previously found in Ref. [lij for 
the case when the hard scale is imposed by the mass of a heavy flavor. 



VII. NUMERICAL RESULTS 

We now turn to a discussion of the numerical results for the most important observables. 
First of all, we are interested in the di-lepton (a la Drell-Yan pair) production channel as 
the simplest one. Although quark production channel could also be of interest, this case will 
be considered elsewhere. 

In Fig. □ (for RHIC energy = 500 GeV) and Fig. |2] (for LHC energy = 14 
TeV) we present the single diffractive cross sections for Z°, 7* (diffractive DY) and 
bosons production, differential in the di-lepton mass squared da s d/dM 2 (left panels) and 
its longitudinal momentum fraction, da s d/dxi (right panels). These plots do not reflect 
particular detector constraints - a thorough analysis including detector acceptances and 
cuts has to be done separately. The M 2 distributions here are integrated over the ad hoc 
interval of fractional boson momentum 0.3 < X\ < 1, corresponding to the forward rapidity 
region (at not extremely large masses). Then the mass distribution is integrated over the 
potentially interesting invariant mass interval 5 < M 2 < 10 5 GeV 2 , and can be easily 
converted into (pseudo) rapidity ones widely used in experimental studies, if necessary. 

The M 2 distributions of the Z° and W bosons clearly demonstrate their resonant behav- 
ior, and in the resonant region significantly exceed the corresponding diffractive Drell-Yan 
component; only for very low masses the 7* contribution becomes important (left panels). 
For X\ distribution, when integrated over low mass and resonant regions, diffractive W + 
and 7* components become comparable to each other, both in shapes and values, whereas 
the W~ and, especially, Z-boson production cross section are noticeably lower (right pan- 
els). Quite naturally, the W~ cross section is (in analogy with the well-known inclusive 
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FIG. 1: Diffractive gauge boson production cross section as function of di-lepton invariant mass 
squared M 2 (left panel) and boson fractional light-cone momentum x\ (right panel) in pp collisions 
at the RHIC energy yfs = 500 GeV. Solid, long-dashed, dashed and dotted curves correspond to 
Z, 7*, W + and W~ bosons, respectively. CTEQ10 PDF parametrization j3l| is used here. 



production) smaller than the W + one due to differences in valence u- and <i-quark densities 
(dominating over sea quarks at large x q ) in the proton, the bosons couple to. So the precise 
measurement of differences in forward diffractive W + and W~ rates would allow to constrain 
quark content of the proton at large x q = X\ja. In Fig. [1] and [2], and in all calculation below 
we have used the most recent CTEQ10 valence/sea quark PDFs parametrization 34j, if not 
declared otherwise. 

In Fig. |3] we show the ratio of the longitudinal (L) to transverse (T) gauge boson po- 
larization contributions to the diffractive production cross section. This ratio is presented 
differentially as function of lepton-pair invariant mass squared {p^ A j dM 2 ) / (aj d / dM 2 ) (left 
panel) and gauge boson fractional momentum [a^J dx\) / [o^J dx-y) (right panel). We see 
that the diffractive gauge bosons production process is always dominated by radiation of 
transversely polarized lepton pairs. The ratio <Jl/&t only slightly depends on M 2 and even 
less on pp cm. energy yfs, so it can be considered as energy independent (which, in fact, can 
be already seen from approximate formulae (14 .9p and f)4.10p ). The longitudinal bosons po- 
larization roughly amounts to 10 % at X\ ~ 0.5 and then steeply falls down at large X\ — > 1 




FIG. 2: The same as in Fig. [TJ but for the LHC energy yfs = 14 TeV. 
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FIG. 3: The ratio of the cross sections of longitudinally (L) to transversely (T) polarized gauge 
bosons, as function of the di-lepton invariant mass squared M 2 (left panel), and boson fractional 
light-cone momentum x\ (right panel). Solid, long-dashed, dashed and dotted curves correspond 
to Z, 7*, W + and W~ bosons, respectively. CTEQ10 PDF parametrization 13411 is used here. 
The ratio depends on pp cm. energy only slightly, by a few percents (cf. Ref. [16() over a vast 
multi-TeV interval, so we neglect it here. 



asymptotically approaching relativistic (massless) bosons case given, due to the gauge in- 
variance, by the transverse polarisation only. Such a behavior turns out to be the same as 
in the inclusive Drell-Yan process [15}. At smaller x\ < 0.6 this ratio becomes the same 
for different bosons, whereas at large x\—±\ the relative contribution of the longitudinally 
polarized photon dominates in corresponding ratios for other bosons. 




FIG. 4: The di-lepton transverse momentum qj_ distribution of the doubly-differential diffractive 
cross section at the LHC energy = 14 TeV at fixed di-lepton invariant mass is shown in the left 
panel. The longitudinal-to-transverse gauge bosons polarisations ratio as a function of the di-lepton 
q± is shown in the right panel. In both panels, the invariant mass is fixed as M = Mz in the Z°, 7* 
production case and as M = M\y in the production case. CTEQ10 PDF parametrization 34] 
is used here. 



From the phenomenological point of view, the distribution of the forward diffractive 
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cross section in the di-lepton transverse momentum q± could also be of major importance 4 . 
In Fig. H] (left panel) we show the di-lepton transverse momentum q± distribution of the 
doubly-differential diffractive cross section at the LHC energy yfs = 14 TeV at the di-lepton 
invariant mass, fixed at a corresponding resonance value - the Z or W mass. The shapes 
turned out to be smooth and the same for different gauge bosons, and are different mostly 
in normalisation. In Fig. [4] (right panel) we show the q± dependence of the o L jo T ratio in 
the resonances. We notice that the ratio does not strongly vary for different bosons. It is 
peaked at about the half of the resonance mass, and uniformly decreases to smaller /larger 
q± values. 
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FIG. 5: Charge asymmetry in the single diffractive W + and W~ cross sections as a function of 
M 2 , at fixed x\ = 0.5 (left panel), and x\, at fixed M 2 = M 2 ^ (right panel). Solid lines correspond 
to the LHC energy yfs = 14 TeV, dished lines - to the RHIC energy y/s = 500 GeV. 

As one of the important observables, sensitive to the difference between u- and <i-quark 
PDFs, the W ± charge asymmetry Aw is shown in Fig. [5] differentially as a function of the 
di-lepton invariant mass squared M 2 and integrated over 0.3 < x% < 1.0 interval (left panel) 



A W {M 2 



daZ + /dM 2 -daZ~/dM 2 



daZ + /dM 2 + daf d -/dM 2 



(7.1) 



and as a function of the boson momentum fraction x\ and integrated over 5 < M 2 < 10 5 
GeV 2 interval (right panel) 



Aw(xx) 



jdx\ 



daf d /dxi 



daZ + /dx 1 + daZ'/d Xl 



(7.2) 



The ratio turns out to be independent on both the hard scale M 2 and the cm. energy 
\fs. One concludes that, due to different x-shapes of valence u, d quark PDFs, at smaller 
X \ < 0.9 the diffractive W + bosons' rate dominates over W~ one. However, at large X\ — > 1 
the W~ boson cross section becomes increasingly important and strongly dominates over 
the W + one. 



4 Authors are indebted to Torbjorn Sjostrand for pointing out this point. 
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FIG. 6: The diffractive-to-inclusive ratio of the gauge bosons production cross sections in pp 
collisions derived in Eq. (|5.6p as a function of the cm. energy y/s (left panel) and the di-lepton 
invariant mass M 2 (right panel). It does not depend on the type of the gauge boson and quark 
PDFs. 



An important feature of the diffractive-to-inclusive Abelian radiation cross sections ratio 



R(M 2 , Xl ) 



da s< il dx\dM 2 
doindl 'dx\dM 2 



(7.3) 



which makes these predictions different from ones obtained in traditional diffractive QCD 



factorisation-based approaches (see e.g. Refs. ll2[), is their unusual energy and scale 
dependence demonstrated in Fig. EJ Notice that we stick to the case of small boson transverse 
momenta, q± <C M, where the main bulk of diffractive signals comes from. The analytic 
formula for this ratio was derived above and is shown by Eq. ( 15. 6ft . which deminstrates that 
the ratio is independent of the type of the gauge boson, its polarisation, or quark PDFs. In 
this respect, it is the most convenient and model independent observable, which is sensitive 
only to the structure of the universal elastic dipole amplitude (or the dipole cross section), 
and can be used as an important probe for the QCD diffractive mechanism for forward 
diffractive reactions, essentially driven by the soft interaction dynamics. We see from Fig. [6] 
that the <y S( i/<J inc i ratio decreases with energy, but increases with the hard scale, thus behaves 
opposite to what is expected in the diffractive factorisation-based approaches. Therefore, 
measurements of the single diffractive gauge boson production cross section, at least, at two 
different energies would provide important information about the interplay between soft 
and hard interactions in QCD, and its role in formation of diffractive excitations and color 
screening effects. 

Finally, in Fig. [7] we present the diffractive-to-inclusive cross sections ratio as a function 
of the boson fractional momentum aq at RHIC (y/s = 500 GeV) and at LHC (y/s = 14 TeV) 
energies (left panel). In the right panel, we compare this ratio calculated at the Tevatron 
energy (y/s = 1.96 TeV) with the recent measurements of diffractive W and Z production 
performed by the CDF collaboration The data show the aq-integrated ratio of diffractive 
to inclusive cross sections. Due to the weak oq-dependence of this ratio (left panel) with a 
good accuracy the integrated values are numerically close to the ratio of differential cross 
sections given by Eq. (17.31) . We see that the results of our calculations with Eq. (15. 6p at 
Xi = 0, 5 agree well with the data. This agreement is another confirmation of correctness of 
the absorption effects included into the parametrization of the dipole cross section (13.121) . 
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FIG. 7: The diffractive-to-inclusive ratio as function of the boson fractional momentum x\. 



VIII. LINK TO THE REGGE PHENOMENOLOGY AND DATA 

The process pp — > Xp at large Feynman xp — > 1 of the recoil proton, or small 



£ = 1 - x F 



x 



(8.1) 



is described by triple Regge graphs as is illustrated in Fig. [HJ Here we also included pro- 



G* J 





FIG. 8: Triple- Regge description of the process pp — > Xp, where the diffractively produced state 
X contains a gauge boson decaying to a lepton pair (see text). 



duction radiation of a gauge boson. The incoming proton can diffractively dissociate either 
due to excitation of the quark skeleton, what corresponds to the triple- Regge graph IPIPM, 
or via gluon radiation corresponding to the triple-Pomeron term, 1P1P1P. So far our calcu- 
lations did not specify which of these terms is related to radiation of the gauge boson. At 
first glance we considered only the valence quark Fock component of the proton. In such 
a case only the IPIPM graph would contribute, what corresponds to the configuration in 
the projectile proton, when the main momentum fraction is carried by two quarks, while 



the third one has a tiny fractional momentum x 



l q,T 



/Mx, where m 2 



m; 



+ p ± is the 



transverse quark mass squared. The gauge boson can be radiated either from the quark 
with minimal fractional momentum (the first term in the r.h.s. of Fig. |H]), or from one of 
the two high fractional momentum quarks (the second term in Fig. [S]). 

The triple-Pomeron term corresponds to diffractive excitation via gluon radiation. This 
process corresponds to the projectile configuration in which the fractional momenta of the 
valence quarks are of the same order, while a gluon (or several of them) carry a tiny fraction. 
In this case the gauge boson can be radiated either by a sea quark (the first term in Fig. [8]), 
or by one of the valence quarks (the second term in Fig. [8]) . 
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Gluons in the light-cone wave function of the proton should also be considered as specta- 
tor partons, and the large distance Rij in Eq. ( 14. ip in this case is the quark-gluon separation. 
In fact, our calculation did include such configurations. Indeed, data on diffraction show 
that diffractive gluon radiation is quite weak (well known smallness of the triple-Pomeron 
coupling), and this can be explained assuming that gluons in the proton are located within 



small "spots" around the valence quarks with radius r ~ 0.3 fm [13|, |37H39| . Therefore, the 
large distance between one valence quark and a satellite-gluon of another quark is approx- 
imately equal (with 10% accuracy) to the quark-quark separation. Since a valence quark 
together with co- moving gluons is a color triplet, in our calculations the interaction am- 
plitude of such an effective ( "constituent" ) quark with the target is a coherent sum of the 
quark-target and gluon-target interaction amplitudes. 

Thus, our calculations effectively cover the gluon radiation, so the triple-Pomeron term 
is included. This is important because this term dominates the diffractive cross section 40 . 
So we can compare with available data from the CDF experiment on W and Z diffractive 
production depicted in Fig. 
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FIG. 9: The diffractive-to-inclusive ratio as function of the invariant mass squared of the produced 
dilepton. The CDF data for W and Z production were taken at the Tevatron energy (y/s = 1.96 
TeV). The first CDF data point corresponds to the W production, M 2 = Myy, the second - to the 
Z production, M 2 = Mf . 



However, in order to compare our results with CDF data, we have to introduce in our 
calculations the accessible experimental cuts, namely, 0.03 < £ < 0.1 [9j. Since our diffrac- 
tive cross section formulae are differential in M 2 , not in M|-, and experimental cuts on 
^/-rapidity distribution of a produced gauge boson are unavailable at the moment, a direct 
implementation of the £ cuts into our formalism cannot be performed immediately. 

As a way out of this problem, at small £ — > one can instead write the single diffractive 
cross section in the phenomenological triple-Regge form [IF 



d£ dp\ 



s ~^J2 6 



I + G 31P (0) e _ BM 



(8.2) 



where s x = 1 GeV 2 ; Bp^i = FP PFi - 2a'jp ln£; i = JP,M; and a'p « 0.25GeV" 2 is the 
slopes of the Pomeron trajectory. Then an effect of the experimental cuts on £ in the 
phenomenological cross section (18. 2p and in our diffractive cross section calculated above 
( 13 .6p should roughly be the same. 
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Since the data show no substantial rise of the diffractive cross section with energy 4l|, |42 



which is apparently caused by strong absorptive corrections, we incorporate this fact fixing 
the effective Pomeron intercept at ajp(O) = 1. This also allows us to use the results of 
the comprehensive triple- Regge analysis of data performed in Ref. (4o|, which led to the 
following values of the parameters: G3jp(0) = Gjpph(0) = 3.2 mb/GeV 2 ; R\jp = 4.2 GeV~ 2 ; 
r wpm = 1-7 GeV -2 . 

Now we are in a position to evaluate the suppression factor 5 caused by experimental cut 
on £, by taking the ratio 

J dp\ J ° 3 d£ da I dv\di 



J dp\ L max d£ da I dp\d^ 



(8.3) 



Here ^ min = Mj [min /s, where M x ,min — M z is the minimal produced diffractive mass 
containing a heavy gauge boson. The value of 5 in Ref. (18.31) is essentially determined by the 
experimental cuts on £ and is not sensitive to the upper limit £ ma z in denominator, so we 
fix it at a realistic value 5 : £ 

max ~ 0.3. Then Ecj. ( I8.3P leads to 5 — 0.2, the factor reducing 
the diffractive gauge bosons production cross section calculated above. Our result plotted 
in Fig. [9] demonstrates a good agreement with the CDF data on single diffractive W and Z 
production [§]. 



IX. CONCLUSIONS 



The diffractive radiation of Abelian fields, 7, Z°, W , expose unusual features, which 
make it very different from diffraction in DIS, and lead to a dramatic breakdown of QCD 
factorisation in diffraction. 

The first, rather obvious source for violation of diffractive factorisation is related to 
absorptive corrections (called sometimes survival probability of large rapidity gaps). The 
absorptive corrections affect differently the diagonal and off-diagonal terms in the hadronic 
current [35|, leading to an unavoidable breakdown of QCD factorisation in processes with 
off-diagonal contributions only. Namely, the absorptive corrections suppress the off-diagonal 
diffraction much stronger than the diagonal channels. In the diffractive Abelian radiation 
in hadron-hadron collisions a new state, i.e. the gauge boson decaying into the heavy lepton 
pair, is produced, hence, the whole process is of entirely off-diagonal nature, whereas in the 
diffractive DIS contains both diagonal and off-diagonal contributions Q|. This is the first 
reason why QCD factorisation is broken in the diffractive gauge bosons production processes. 

The second, more sophisticated reason to contradict diffractive factorisation is specific for 
Abelian radiation, namely, a quark cannot radiate in the forward direction (zero momentum 
transfer), where diffractive cross sections usually have a maximum. Forward diffraction 
becomes possible due to intrinsic transverse motion of quarks inside the proton. 

Third, the mechanism of Abelian radiation in the forward direction in pp collisions is 
related to participation of the spectator partons in the proton. Namely, the perturbative 



5 The estimate for ^ max ~ 0.3 corresponds to the limiting case when one of the constituent quarks in a 
target (anti)proton looses almost all its energy into a hard radiation of a gluon in the t-channel. The 
second and all subsequent i-channel gluon exchanges collectively screen the color charge taken away from 
the target by the first gluon and transfer much smaller fraction of initial target momentum to projectile 
quarks as has been recently advocated in Refs. 



43, 



44|. 
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QCD interaction of a projectile quark is responsible for the hard process of a heavy boson 
radiation, while a soft interaction with the projectile spectator partons provides color neu- 



tralization |43|, |44j, which is required for a diffractive (Pomeron exchange) process. Such 
an interplay of hard and soft dynamics is also specific for the process under consideration, 
which makes it different from the diffractive DIS, dominated exclusively by soft interactions, 
and which also results in breakdown of diffractive factorisation. 

The diffractive (Ingelman-Schlein) QCD factorisation breaking manifests itself in specific 
features of diffractive observables like a significant damping of the single diffractive gauge 
bosons production cross section at high y/s compared to the inclusive production case. This 
is rather unusual, since a diffractive cross section, which is proportional to the dipole cross 
section squared, could be expected to rise with energy steeper than the total inclusive cross 
section, like it occurs in the diffractive DIS process. At the same time, the ratio of the single 
diffractive to inclusive production cross sections rises with the hard scale, M 2 . This is also 
in variance with diffraction in DIS associated with the soft interactions. 

In this paper, we have presented the differential distributions (in transverse momentum, 
invariant mass and longitudinal momentum fraction) of the diffractive 7*, Z° and W ± bosons 
production at RHIC (500 GeV) and LHC (14 TeV) energies, as well as the ratio of the boson 
longitudinal to transverse polarisation contributions. We have also calculated the charge 
W ± asymmetry, relevant for upcoming measurements at the LHC. The ratio diffractive to 
inclusive gauge bosons production cross sections does not depend on a particular type of 
the gauge boson, its polarisation state and quark PDFs, and depends only on properties of 
the universal dipole cross section and sensitive to the saturation scale at small x. Finally, 
our prediction for this ratio is numerically consistent with the one measured for diffractive 
W and Z production at the Tevatron. 

The theoretical uncertainties of our calculations come mainly from poorly known quark 
PDFs (or the proton structure function F2(x q ) in diffractive DY) at large quark fractions 
x q — y 1. This issue has been discussed in our previous analysis of the DDY process in 
Ref. [l6j, where a strong sensitivity of the x\ dependence to a particular F2 parametrization 
has been pointed out. The same situation extends to the more general case of diffractive 
Abelian radiation considered in this paper. This tells us again that measurements of 
forward diffractive gauge bosons production would be extremely important or even crucial 
for settling further more stringent constraints on the quark content of the proton. 
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